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Abstract 

In this paper we study two multicritical correlation kernels and prove that they 
converge to the Pearcey kernel in a certain double scaling limit. The first kernel 
appears in a model of non-intersecting Brownian motions at a tacnode. The second 
, arises as a triple scaling limit of the eigenvalue correlation kernel in the Hermitian 

Qh two-matrix model with quartic/quadratic potentials. The two kernels are different 

but can be expressed in terms of the same tacnode Riemann-Hilbert problem. The 
proof is based on a steepest descent analysis of this Riemann-Hilbert problem. A 
special feature in the analysis is the introduction of an explicit meromorphic function 
on a Riemann surface with specified sheet structure. 
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1 Introduction 

Point processes with determinantal correlation kernels have attracted a lot of interest 
over the past few decades due to the rich mathematical structures behind them and 
their frequent occurrences in various random models including random matrix theory, 
random growth and tiling problems, etc. [34|, 139]. A fundamental issue in the study of 
determinantal point processes is to establish microscopic limits of the correlation kernels, 
which often leads to universal results. Besides the well-known canonical kernels like the 
sine kernels for bulk universality, the Airy/Bessel kernels for soft /hard edge universality 
and so on, some new kernels were found describing critical behavior of certain random 
models, and they are believed to be universal as well. This paper deals with two such 
models in which a similar phenomenon occurs. We intend to retrieve a canonical process 
(the Pearcey process here) from the critical phenomenon. Physically, this will lead 
to descriptions of phase transitions among different processes. For similar transitions 
between canonical processes, we refer to the thesis of Deschout [25J and the recent papers 



1.1 Non-intersecting Brownian motions at a tacnode 

As a first model we consider n one-dimensional non-intersecting Brownian motions with 
two starting points at time r = and two ending points at time r = 1. The transition 
probability density of the Brownian motions is given by 



mum- 




(1.1) 
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Figure 1: Non-intersecting Brownian motions with two starting points ±a and two 
ending positions ±/3 in case of (a) large, (b) small, and (c) critical separation between 
the endpoints. Here the horizontal axis denotes the time, r G [0, 1], and for each fixed r 
the positions of the n non-intersecting Brownian motions at time r are denoted on the 
vertical line through r. Note that for n — > oo the positions of the Brownian motions fill 
a prescribed region in the time-space plane, which is bounded by the boldface lines in 
the figures. Here we have chosen n = 20, T = 1 in each of the figures, and (a) a = 1, 
(3 = 0.7, (b) a = 0.4, /3 = 0.3, and (c) a = 1, /3 = 0.5, in the cases of large, small and 
critical separation, respectively. 

where we interpret T > as a temperature variable. As the number of paths tends to 
infinity, these paths will fill out a certain domain in the time-space plane. Depending on 
the locations of the starting and ending points, the temperature T, and the fractions of 
paths connecting the topmost and bottommost starting and ending points, we distinguish 
three cases, namely, large, small, and critical separation. This is illustrated in Figure Q] 
which is taken from [IT], see also [23] . 

For each fixed r £ (0, 1), the positions x\, . . . , x n of the Brownian paths form a deter- 
minantal point process. It is well-known that the scaling limits of the correlation kernel 
are given by the sine kernel in the interior of the domain, the Airy kernel at a typical 
point of the boundary, and the Pearcey kernel |13t 114] at a cusp |12[ 140]. In the critical 
situation the Brownian paths fill out two touching ellipses; see Figure HJc). The point 
where the ellipses touch is called the tacnode and the local particle correlations around 
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this point are described by the tacnode kernel. This tacnode process was recently studied 
by different groups of authors using different techniques. Adler-Ferrari-Van Moerbeke 
[3] resolved the tacnode problem for non-intersecting random walks (discrete space and 
continuous time). Johansson [35] gave an integral representation of the extended tacn- 
ode kernel in the continuous time-space setting; see also Ferrari and Veto [31] for the 
asymmetric tacnode process. Delvaux-Kuijlaars-Zhang [23] found an expression of the 
tacnode kernel in terms of a new 4x4 Riemann-Hilbert (RH) problem. Furthermore, the 
tacnode process also appears in a random tiling model called the double Aztec diamond 

The critical separation case, see Figure QJc), can be seen as a limit of the small 
separation case, in which the two groups of particles merge at a certain time and separate 
again later, see Figure [T^b). In this case the particles fill out a certain shape in the time- 
space plane with two cusp points. At these cusp points the limiting particle density 
vanishes with an exponent 1/3 which indicates that the local particle correlations are 
governed by a Pearcey kernel. When we increase the distance between both starting 
points and both ending points keeping the temperature T fixed, or fix the starting and 
ending points while lowering the temperature T, the cusp singularities approach each 
other to form a tacnode. Hence we expect to retrieve the Pearcey kernel from the tacnode 
kernel in a certain limit. One aim of this paper is to describe such a transition. 

1.2 Hermitian two-matrix model with quartic/quadratic potentials 

A similar phenomenon occurs in the Hermitian two-matrix model, which is a probability 
measure 

— exp (-nTr(y(Mi) + W(M 2 ) - tMM)) dMi dM 2 , (1.2) 

defined on the space of pairs (Mi, M%) of Ti x n Hermitian matrices. Here, Z n is a nor- 
malization constant, r > is the coupling constant, and V and W are two polynomials 
which we take as 

V(x) = ^, and W(y) = ^ + a^, (1.3) 

for a parameter a£l. 

The two-matrix model can be integrated in terms of biorthogonal polynomials, see 
[291 El [TU] . The correlation functions for the eigenvalues of M\ and M 2 have a deter- 
minantal structure and admit expressions in terms of these polynomials [3Q|. [36] . Here 
we are only interested in the eigenvalues of M\ when averaged over M 2 , which form a 
determinantal point process. 

As the parameters a and r vary, we encounter phase transitions. In particular, the 
point (—1, 1) in the ar-plane corresponds to a multicritical case. A new kernel, referred 
to as the critical kernel in this paper, was recently established by Duits and Geudens 
[26j to describe the local eigenvalue correlations near this point. This critical kernel is 
expressed in terms of a RH problem (see RH problem 11.11 below) that is a generalization 
of the one used in the tacnode model [24] . Nevertheless, the critical kernel is genuinely 
different from the tacnode kernel. Although at the first glance not so clear as in the 
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non-intersecting Brownian motions model, an analysis of the ar-phase diagram reveals 
that this critical kernel can be understood as a limiting case of a Pearcey kernel or a 
Painleve II kernel 1 15 j : see Section [3.11 for details. In |26j it was already shown that 
the critical kernel reduces to the Painleve II kernel by taking a double scaling limit. The 
other aim of this paper is then to retrieve the Pearcey kernel from the critical kernel, 
hence, complementing the results in [26j. 

1.3 The tacnode Riemann-Hilbert problem 

The essential feature allowing us to handle the tacnode kernel and the critical kernel 
simultaneously is that both kernels have explicit representations in terms of the following 
tacnode RH problem [221126]. 

RH problem 1.1 (Tacnode RH problem). Fix parameters ri,r 2 , si, s 2 , andt. We look 
for a 4 x 4 matrix-valued function M(z) satisfying 

(1) M is analytic for z £ C \ Em- 

(2) For z £ Tfc, the limiting values 



exist, where the -\--side and —-side of r& are the sides which lie on the left and 
right ofT^, respectively, when traversing according to its orientation. These 
limiting values satisfy the jump relation 



M+(z) 



lim 

f on +-side of I\ 



M(C) 



Af_ (z) 



lim 

(on — side of Tj. 



M(C) 



M+(z) = M.{z)J k {z) 



fc = 0,...,9. 



(1.4) 



(3) As z — > oo with z £ C \ Y<m we have 



M{z) = (I + Oiz- 1 )) B{z)A 




where 



Mz) = 2 -r x z^ 2 + 2 Sl z 1 ' 2 , ^ 2 {z) = \r 2 (-zfl 2 + 2s 2 {-z)V 2 



(1.6) 
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Figure 2: The jump contour Em in the complex z-plane and the constant jump matrices 
Jfc on each of the rays T^, k = 0, . . . , 9. 



(4) M(z) is bounded near z = 0. 

The contour £ m is shown in Figure and consists of 10 rays emanating from the 
origin. The function M(z) makes constant jumps on each of the rays IV These rays 
are determined by two angles tp\ and (p2 satisfying < cpi < (f2 < vr/2. The half-lines 
Tfz, k = 0, . . . , 9, are defined by 

T = [0,oo), ri = [0,e^oo), T 2 = [0,e^ 2 oo), T 3 = [O^^cjo), T 4 = [0, e^-^oo), 
and 

T5+k = -Tk, fc = 0, ...,4. 

^4ZZ rays are oriented towards infinity. 

The fractional powers are defined with respect to the principal branch. Hence for 
example, z i-> (— z) 3 / 2 is analytic inC \ [0,oo) and takes positive values on the negative 
part of the real line. 

By [211 [26], this RH problem has a unique solution for r\ = r 2 > 0, si = s 2 G M, 
and t £ i, Moreover the tacnode RH problem has a remarkable connection with the 
Hastings-McLeod solution of the Painleve II equation [33] . 

In Sections 12.21 and 13.21 below, we will see that the tacnode RH problem appears in 
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expressions of the tacnode and critical kernel for the choices of parameters 

'n =7-2 = 1, 

< si = s 2 = s e R, (1.9) 
k t G R. 

We denote the associated solution as M(z;s,t). 
1.4 Outline of the paper 

The rest of this paper is organized as follows. Our results for the non-intersecting 
Brownian motions model and the Hermitian two-matrix model are stated in Sections [2] 
and [3l respectively. In both sections, we first discuss the model from the angle of the 
phase diagram, which is particularly helpful in understanding the critical phenomena and 
phase transitions that occur. Next we give explicit expressions of the tacnode and critical 
kernel in terms of the unique solution to RH problem 11.11 with parameters (|1.9p . Our 
main results are Theorems 12.31 and 13.31 which state that the tacnode and critical kernel 
converge to the Pearcey kernel in a certain double scaling limit. The proofs are given 
in Section [6j They are based on the Deift-Zhou steepest-descent analysis [HI QjJl [20] of 
RH problem ll.il which will be performed in Section [SJ A special feature in the analysis 
is the introduction of an explicit meromorphic A-function on a Riemann surface with 
specified sheet structure, which will be the topic of Section [H We emphasize that our 
approach is based on the steepest-descent analysis of a larger size RH problem as in [25] 
and is different from [2j |6l [7] . 



2 From the tacnode kernel to the Pearcey kernel 
2.1 The phase diagram 

We consider a symmetric version of a model of n 1-dimensional non-intersecting Brow- 
nian motions as in |1T|. 123], The Brownian particles start at points rtor, a > 0, at time 
r = and end at points ±/3, /3 > 0, at time r = 1. More precisely we assume that 
n is even and that n/2 of the particles move from the topmost starting point a to the 
topmost ending point (3, while the other half of the particles connect the point —a to 
the point — /3. 

For convenience we fix the starting and ending points such that 

2a(3 = 1. (2.1) 

According to [23], the three cases of large, small, and critical separation of the starting 
and ending points then correspond to T < 1, T > 1, and T = 1, respectively; see Figure 
[1] again. Moreover, in case of critical separation, the coordinate of the tacnode is given 
by (r cri t,0), where 

r crit = ^ (2.2) 
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Figure 3: Phase diagram for n non-intersecting Brownian motions in the symmetric 
setting, n/2 particles connect a and f3 while the other particles move from —a to —(3. 
In this picture a = (3 = l/v2- The boldface curves indicate phase transitions and have 
equations T = 1 and T = (r 2 — r + l/2)/(r(l - r)). 

is the critical time. 

The assumption (|2.1|) allows us to summarize the situation in a clear phase diagram, 
shown in Figure [3] in the case a = (3 = 1/V2. There are two curves in the TT-plane. 
One is the straight line T = 1. The other one is given by the algebraic equation 

T{t) = a 2 — + /3 2 -^. (2.3) 
r 1 — r 

These two curves divide the phase plane into three regions, which we denote by Case I, 
Case II, and Case III. These three cases correspond to three generic situations and the 
two curves indicate phase transitions. Case I corresponds to temperature T < 1, i.e., the 
case of large separation. In this case, the Brownian paths remain in two separate groups, 
and the limiting hull in the rx-plane consists of two disjoint ellipses, see FigureQJa). The 
particles are asymptotically distributed on two intervals according to semicircle laws. In 
Cases II and III we have T > 1 with the paths being distributed on a single interval 
in Case III and on two intervals in Case II. In both cases the densities of the limiting 
distribution for the positions of the paths at time r are different from semicircle laws. 
The curve T = 1, which corresponds to critical separation in Figure [U separates Case 
I from Case II. The transition occurring at this curve was studied in [23 and is related 
to the homogeneous Painleve II equation. This phase transition is not visible in the 
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correlation kernel, but manifests itself in critical limits of the recurrence coefficients of 
the associated multiple orthogonal polynomials. The curve (|2.3[) separates Case II from 
Case III. This transition corresponds to the cusps in the limiting hull in Figure QJb), 
which indicates local behavior described by Pearcey kernels. 

The two critical curves touch at the multicritical point T = 1, r = r cr it, which 
corresponds to the tacnode point in the critical separation. At this point the limiting 
particle density consists of two semicircles touching at the origin. In [24J the local 
particle correlations around the origin were studied for fixed temperature T = 1 and 
time r = r cr it, i.e. exactly at the multicritical point (l,r cr i t ), and varying starting and 
ending points. Under these assumptions the local scaling limit is expressed in terms 
of RH problem 11.11 with the parameter t = 0. In this paper, however, we find it more 
convenient to fix the starting and ending points a and j3 according to (|2.ip and scale 
(T, r) around the multicritical point as follows: 

r = r crit + Kn- 1 / 3 = + Kn- l '\ 

a + P (2.4) 

T = 1 + Ln~ 2/3 , 

where K and L are arbitrary real constants; see also [22] for a similar scaling. Then, 
as n — > oo, the correlation kernel of non-intersection Brownian motions K n (x,y) at the 
tacnode will converge to the tacnode kernel K t£ic introduced in the next section, after 
proper scalings of the space variables. 



2.2 The tacnode kernel 

Definition 2.1. For u,d£1, we define the tacnode kernel K t&c by 



K tac (u,v;s,t) 



1 



2m{u — v) 



(0 1 1) M _1 (t;;s,t)M(u;s,t) 



i 





(2.5) 



where M{z; s, t) denotes the analytic continuation of the restriction of M(z; s, t) to the 
sector around the positive imaginary axis, bounded by the rays T2 and Recall that 
M(z; s,t) is the unique solution to RH problem 11.11 with parameters (|1 .9[) . 

We can rewrite the kernel in terms of the limiting values M + (n; s, t) and M + (v; s, t) 
of M on the real line, by using the jump relations in RH problem ll.il For example, for 
u, v > 0, we have 

A\ 

K tac (u,v;s,t) = — -1 -(-1 1 0) M^(v; S ,t)M + (u; S ,t) ° 

lm{u — v) 1 

with a different expression in case u and/or v are negative. 

We then have the following results concerning the correlation kernel of non-intersecting 
Brownian motions at a tacnode. 



(2.6) 



9 



Theorem 2.2. Consider n nonintersecting Brownian motions on R with two starting 
points ±a, a > and too ending points ±/3, /3 > 0. The transition probability density is 
given by (jl.ip . Suppose half of the paths start in a at time r = and end in f3 at time 
r = 1, while the other half of the paths connect —a and —(3. If we fix a and f3 such that 
2af3 = 1 and make a double scaling limit of the temperature T and the time r according 
to (|2.4p . the correlation kernel K n for the positions of the particles near the critical time 
satisfies 

lim — TTT^Kr, f 



= K tac (u, v- 2- 5 / 3 ((a + f3) 4 K 2 - L) , -2~ 1 / 3 (a + (3) 2 K), (2.7) 
uniformly for u,v in compact subsets ofM, where c = 2 1 / 3 (a + /3) and K tac is given by 

dm 



Proof. A version of this theorem with different scaling assumptions was proved in [24] . 
The proof of this theorem follows the same lines. The main difference is in the construc- 
tion of the local parametrix around the origin, which is now built out of the solution to 
RH problem 11.11 for nonzero i, whereas in |24j the version for t = was used. See also 
22 1 for a similar construction in a more involved model. □ 



2.3 Double scaling limit: from the tacnode kernel to the Pearcey kernel 

The Pearcey kernel K Pe is defined as (see |13t 114]) 

k Pc (x y p) = E^KM -p'( x )i'(y) +p"( x )i(y) - pp{x)g{y) , 2 g s 

x-y 

with 

p(x) = — [°° e~\ sA -2 s2+isx ds and q(y) = — [ e^ t4+ 2 t2+ity dt, (2.9) 
2vr y.oo 27r J s 

where pel. The contour S consists of the four rays argi = 7r/4, 3tt/4, 57r/4, 77r/4, 
where the first and the third ray are oriented from infinity to zero while the second and 
the last ray are oriented outwards. The functions (|2.9|) are called Pearcey integrals |38j 
and are solutions of the third order differential equations p"'{x) = xp(x) + pp{x) and 

Q"'(y) = -M{y) + PQ(y), respectively. 

The Pearcey kernel appears in several random models including the closing of gap 
in the Gaussian random matrix model with external source [13} \12\ 00] , non- intersecting 
Brownian motions at cusps [5] (see also interpretations in this paper), a combinatorial 
model on random partitions [37] . etc. We mention that, besides the expression (|2.8p . the 
Pearcey kernel admits a double integral representation as well as a RH characterization. 
For the latter one, see Section 16.11 

This is our first main result. 
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Theorem 2.3. Let K tAC be as defined in (|23|) and K Pe as in (|2~%|) . T/ien 



-oo 



lim -= -K tac -= » • ,| a | 1 + — — = i^ e (y,x;cT) 

-V^|a|V4 \V2\a\V* V2|a|V4' 2 ' 11 V 2|o| 3 / 2 , H 



and 



(2.10) 

lim * ; if tac f J* ; , J* ; ■- — ,- a (l + -^ 7 -) \ = K Fe (x,y;a), (2.11) 

uniformly for x, y in compact sets, where a 6l is a parameter. 

The good scaling of the parameters s,t in the tacnode kernel K tac can be derived 
from the phase diagram in Figure [3l Indeed, the phase diagram suggests that it is 
possible to retrieve the Pearcey kernel from the tacnode kernel by letting r and T move 
away from the multicritical point (r cr it, 1) along the curve (j2.3|) that separates Case II 
from Case III. Recall the double scaling (|2.4|) of r and T involving K and L. If we 
leave the multicritical point to the left/right along the curve, by (|2,4p . this amounts to 
setting L = 2K 2 (a + /3) 4 and =F-fT > 0. This implies that the tacnode kernel appearing 
in Theorem 12.21 should take parameters s,t as 

8 = -a 2 /2 and t = -a, (2.12) 

with =Fa = =p2 _1 / 3 (a+/3) 2 i^ > 0. A further consideration noting that the Pearcey kernel 
is a one parameter family invokes us to take 

f = -a (l + ^) , o-GM, (2.13) 

in our setting. 

Theorem 12.31 will be proved in Section 16.21 In Sections SHS we will mainly focus on 
proving (|2.10p . In the end we will obtain (|2.1ip from (|2.10|) using a symmetry argument. 
We next state the results for the transitions in the Hermitian two-matrix model. 



3 From the critical kernel to the Pearcey kernel 
3.1 The phase diagram 

As aforementioned, the eigenvalues of the matrix Mi in the Hermitian two-matrix model 
(|1.2p form a determinantal process with correlation kernel Kj^\x,y). In this section 
we will again use a phase diagram to get more insight into the nature of the critical 
phenomena and phase transitions involved. 

By the results in [28], it follows that the mean eigenvalue density for Mi has a limit 
as n — > oo, i.e., there exists an absolutely continuous measure hi on R such that 

lim -K[^\x,x) = ^. (3.1) 
n->oo n ax 
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For technical reasons this result was only proved under the assumption that n = 
mod 3. A crucial point in the analysis of [28] is that the limiting measure can be 
characterized by a vector equilibrium problem for vectors of three measures (/^i, /U2, ^3), 
see also [32]. This equilibrium problem involves a constraint 02 acting on the second 
measure fi^i where 02 is a certain measure on the imaginary axis. Moreover, the supports 
of the measures fi±, 02 — /X2, and /i 3 have the following forms 

suppOi) = [-fa,-fa]U [fa, fa], 
supp(<T 2 - M 2 ) = iR \ (-ifa,ifa), 
supp(/i 3 ) = R\(-/3 3 ,/3 3 ), 

for some /3o > /3i > 0, fa, fa > that all depend on the parameters a S R and r > 0. 
Generically, at least one of these numbers is zero, and no two consecutive ones are zero. 
This leads to the following classification of cases as proved in [27\ |2"8] . 

Case I: fa = 0, fa > 0, and fa = 0. In this case there are no gaps in the supports 
of the measures [i\ and // 3 on the real line. The constraint 02 is active along an 
interval [—ifa,ifa] on the imaginary axis. 

Case II: fa > 0, fa > 0, and fa = 0. In Case II there is a gap in the support of fj,%, 
but there is no gap in the support of // 3 , which is again the full real line. The 
constraint is active along an interval along the imaginary axis. 

Case III: fa > 0, fa = 0, and fa > 0. In Case III there is a gap in the supports of \x\ 
and /i 3 , but the constraint on the imaginary axis is not active. 

Case IV: fa = 0, fa > 0, and fa > 0. In this case the measure [i\ is still supported on 
one interval. However there is a gap (— fa, fa) in the support of ,u 3 . As in Case I, 
the constraint 02 is active along an interval [— ifa,ifa] on the imaginary axis. 

The correspondences between the values of (a, r) and these four cases are illustrated 
in the phase diagram shown in Figure HI The different cases are separated by two curves 
with equations 

r = Va + 2, -2 < a < 00, (3.2) 

and 

t = \—, -oo<a<0, (3.3) 
V a 

respectively. On these curves two of the numbers fa, fa, and fa are equal to zero. For 
example, on the curve (|3.3|) that separates Case III from Case IV, we have fa = fa = 0, 
while fa > 0. Finally, note the multi-critical point (a, r) = (—1, 1) in the phase diagram, 
where fa = fa = fa = 0. All four cases come together at this point in the QT-plane. 

As long as we consider points (a, r) that are not on the curves, the local eigenvalue 
correlations are governed by the sine kernel in the bulk of the spectrum and the Airy 
kernel at the edge of the spectrum. Critical phenomena occur at the curves that separate 
the different cases. When we cross the line r = yja + 2 for a > — 1, the support of jjL\ 
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Figure 4: The phase diagram in the ar-plane: the critical curves r = y/a + 2 and 
r = yj— i separate the four cases. The cases are distinguished by the fact whether is 
in the support of the measures fi±, o<i — fj>2, and n$, or not. 



turns from two intervals into one interval. On this critical curve, the intervals meet at 
the origin and there the density vanishes quadratically. This indicates that in a double 
scaling limit the correlation kernel converges to the Painleve II kernel that arises at the 
closing of a gap in unitary ensembles [111 [T5"l [T6] . When crossing the line r = -J— 1/a 
for a < —1, we again have a transition of two intervals merging at the origin. However, 
due to the fact that supp(o"2 — ^2) also closes simultaneously, the vanishing at the origin 
occurs with an exponent 1/3. This indicates that the local correlations are governed by 
the Pearcey kernel. The other transitions, represented by the dashed lines in Figure HI 
do not concern They take place on the non-physical sheets of the spectral curve and, 
therefore, they do not influence the local correlations of the eigenvalues of Mi, which 
are again described by the sine and Airy kernels. 

The limiting process near the multicritical point (a, r) = (—1,1), where there is 
a simultaneous transition in the supports of all three measures fix, — ^2 and /i3 is 
the subject of the recent paper [26]. Here the limiting eigenvalue density vanishes like 
a square root at the interior of the support, which is similar to the situation in non- 
intersecting Brownian motions at a tacnode. By taking triple scaling limits such that a 
and t are made dependent on n as 

(;) = (- 1 ) + „-V.0) +h .-^(->), (3.4) 

for a, b G M, the large n limit of the correlation kernel K^\x, y) is given by the critical 
kernel K CT that we will introduce in the next section. Note that in (|3.4p the vectors 
(2 l) and (—1 2) are respectively tangent and normal to both critical curves in the 
point (a,r) = (—1,1). 
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3.2 The critical kernel 

Definition 3.1. For u,v E K, we define the critical kernel K cv by 



K a '(u,v;s,t) = -(-1 1 0) M(m;s,£) _1 M(w;s,f) 

Z7TZ(U — Uj 



i 



w 



(3.5) 



The following theorem is the main result in [26J. 



Theorem 3.2. i26{ Theorem 2.3] Let be the kernel describing the eigenvalues of 
M\ when averaged over M2 in the Hermitian two-matrix model Q1.2]) with V and W as 
in (jl.3p . Assume a and r depend on n as in (|3.4p . Then for n — > 00 and n = mod 6 ; 



uniformly for u,v in compact subsets o/R ; where K cv is defined in (|3,5p 



3.3 Double scaling limit: from the critical kernel to the Pearcey kernel 

Our second main result is the following. 

Theorem 3.3. Let K CT be as defined in (|3"3|) and K Pe as in ([23]) . T/ien 



lim * ; K cr ( ^ X - , J 1 - ;-±a 2 ,\a\ ( 1-—^- ) ] = if Pe (x, y; a), 
^-oov^| a |i/4 VV^Ial 1 ^' ^2| a |i/4' 2 y 2|a| 3 /Vy V ' y ' ; ' 

uniformly for x, y in compact sets, where a £M is a parameter. 

The scaling of the parameters s, t in the critical kernel K cr is again motivated by the 
phase diagram in Figure [H This phase diagram suggests that it is possible to retrieve the 
Pearcey kernel from the critical kernel by letting t and a move away from the multicritical 
point along the curve r = ^/—l/a, a < — 1. This amounts to setting b = 3a 2 / 5 and 
a < 0. Hence the critical kernel in Theorem 13. 21 1: akes parameters s = —a 2 /2 and t = \a\. 
To obtain the limit in the most general form, we perform the scaling for t as 



t = |o| l 1 -2hFj' aeR - (3 ' 6) 

We also note that if one lets t and a move away from the multicritical point along the 
curve r = \fu + 2, a > — 1, the critical kernel will tend to the Painleve II kernel; see 
[26j for more details. 

Theorem 13.31 will be proved in Section! 
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4 Meromorphic A-functions on a Riemann surface 



In this section, we introduce some auxiliary functions and study their properties. The 
aim is to construct the so-called A-functions, of which the analytic continuation defines 
a meromorphic function on a Riemann surface with specified sheet structure. The A- 
functions have desired behavior around each branch point, and we shall make use of 
them in the normalization step of our steepest descent analysis in Section 15.31 

Throughout this section, unless specified differently, we shall take the principal cut 
for all fractional powers. 

4.1 A four-sheeted Riemann surface and the ^'-functions 

We introduce a four-sheeted Riemann surface TZ with sheets 

Ki = C\ [0, oo), 1Z 3 = C \ ([-ic, ic] U [0, oo)) , 



where 



Tl 2 = C \ (-oo,0], TZi = C\([-ic,ic] U (-00, 0]), 

:!n/ V 3/2 , (4.1) 



16 

and 7 > is a parameter that we will specify later. 

We connect the sheets IZj, j = 1, 2, 3, 4, to each other in the usual crosswise manner 
along the cuts [—ic, ic], [0, 00), and (—00, 0]. More precisely, TZ± is connected to 7^3 along 
the cut [0,oo), IZ2 is connected to IZ4 along the cut (— 00, 0], and IZ3 is connected to 
72-4 via [— ic, ic]. Moreover, the Riemann surface is compactified by adding two points at 
infinity. The first point 00 1 is added to the first and the third sheet, while the second 
point 002 connects the other two sheets. This Riemann surface TZ has genus zero and is 
shown in Figure 

We intend to find functions Aj on these sheets, such that each Xj is analytic on 
IZj and admits an analytic continuation across the cuts. Later we will impose extra 
conditions on the poles and zeros of this function. For the moment, we content ourselves 
constructing an elementary function w that is meromorphic on TZ. 

Consider the algebraic equation 

(7 3 ui 2 — l) 2 ^ 

It turns out this equation defines a meromorphic function on the Riemann surface 7Z. 
This is made precise in the following lemma. 

Lemma 4.1. There exist functions Wj(z), j = 1,2,3,4, that solve the equation (|4.2p 
and satisfy the following conditions. 

(a) Wj(z) is analytic on IZj, j = 1,2,3,4, and 

wi,±(x) = w 3)T (x), x G (0, 00), 

W2,±{x) = ui 4iT (x), x e (— oo,0), 

i03,±(z) = W4 tT (z), z G (— ic, 0) U (0,ic). 
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Figure 5: Riemann surface 1Z. 



Here, we orient (— oo,0) and (0, oo) from left to right, and (— ic, 0) U (0, ic) from 
bottom to top. Hence, the function Uj=i ~^ ^ : ^ z i— >■ Wj(z) has an analytic 
continuation to a meromorphic function w : 1Z — >■ C. T/iis function is a bijection. 



(b) w satisfies the symmetry properties 



Wj (z) = Wj (z), 
wi(-z) = -w 2 (z), 
m{-z) = ~wa{z), 



zeKj, j = 1,2,3,4, 

z € C\ (-oo,0], 

z £ C \ ((-oo,0] U [-ic,ic]) . 



(c) >ls z — > oo within C \ (— oo,0] we /tawe that 

W2 ( Z ) = r V2 + h~ 9/4 z- 1/2 - m7" 15/4 ^ 3/2 



f^7- 9/ V 2 - ^7- 21/4 ^ 5/2 + 0(^ 7/2 ), 



„, 4 (s) = 7 "3/ 2 - I 7 - 9 /S" 1 /2 + 1 7 -15/4 z -3/2 



64 
+ 



128 



7 



-9/2, 



+ ll67- 21/4 ^ 5/2 + 0(- 



-7/2 > 



(d) ^4s z — > we have that 
w 2 {z) 



7 -2^-l/3 + | 7 -l z l/3 _ l z + 28^5/3 



0(z 



7/3 > 



in C \ (— oo, 0], 



w 4 (z) 



2 7 3 z 3 + 9 7 6 z 5 + C>(z 7 ), 



in J U JV, 



-U! 



2 7 -2(_ z) -l/3 _ 2 W7 -l(_^l/3 _ l z _ 28 w 2 7( _ z) 5/3 + 0(^7/3^ m jj 



[o; 2 7 - 2 ^" 1/3 + §W7"^ 1/3 " + ff^V 5/3 + 0(^ 7/3 



in JJJ. 



-ffere, /, II, III and IV stand for the four open quadrants in the complex plane and 
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(e) As z — > ic we have that 

Ws{z) = A 7 "3/2 + g . 3 -7/4 e 3*V4 7 -3/4 (j8 _ -^1/2 + Q {z _ ^ > 

v3 

where we take the branch cut of (z—ic) 1 ^ 2 along (— too, ic) and —n/2 < arg(z-ic) < 

2 ■ 

(f) We have 

w 3 , + ((0, ic)) = (0, ^7~ 3/2 ) , ((0, ic)) = (ioc, ^7" 3/2 ) • 

Note that the complete behavior of Wj, j = 1,2,3,4, around the origin, oo, and ±ic 
can be obtained from the symmetry conditions in (b). 

Proof. We consider the algebraic equation defined by 

{w 2 + 7 3 ) 2 = zw 3 , (4.3) 

where again 7 is a positive parameter as in (|4.2|) . For any solution w{z) of (|4,3|) . the 
function 

w{z) = i {w(i/ z))^ 1 

satisfies equation (|4.2p . Our strategy is then to construct Wj from the solutions of (|4.3p 
via the above transformation. Fortunately, the equation (|4.3[) was well-studied in |26} 
Section 3]. For each z E C, it was shown that the equation (j4.3j) admits four solutions, 
denoted by Wj(z), j = 1, 2, 3, 4, such that 

> 1^2(2)1 > 1^3(^)1 > |W4(2)|. 

These functions can be interpreted as a meromorphic function on another four-sheeted 
Riemann surface 7Z with sheets 

n 1 = C\[-c,c], 7Z 3 = C\(RUiR), 

Tl 2 = C \ ([-c, c] U iR) , £ 4 = C\R, 



where 



These four sheets are connected in the usual crosswise manner and the Riemann surface 
is compactified by adding two points at infinity: one on the first sheet, and another to 
the three remaining sheets. Each Wj{z) is analytic on TZj and the function Uj=i^j ~~ ^ 
C : TZj 3 z 1 — y Wj{z) has an analytic continuation to a meromorphic function w : TZ — > C. 
This fact invokes us to define the following functions Wj(z): 

/n /*(w3(«V*or\ in I U TV, /i^t"/*)) -1 , in/U/y 

U>1 (2 = < 1 Worn = < 1 

\i(w 4 (i/z)) ~\ in JJ U III, ~ \i(w 3 (i/z)) , in II U III, 

, [i{w 2 {i/z))- 1 , in/U/y fi^Ci/z)) -1 , in/U/y 

I %(w\(i/z))~ x , in II \J III, \i(w 2 (i/z)y , in II U III. 
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Then, the lemma follows from |26t Lemmas 3.1-3.3] and straightforward computations. 

□ 



4.2 The A-functions 

With the functions Wj(z) in Lemma |4.1|. we define the A-functions as 

Xj(z) = z 2 {C lWj {zf + C 2 Wj(z) 2 + C 3 + C 4Wj {z)- 2 ) , z e Kj, j = l, 2, 3, 4, 



(4.4) 



where 

d = -3^,33/4 , J. 27/4 Q 7 21/4 , 17 15/4 

Co = -5^9/4 _ 65 3/4 ^ = 1-3/4 , ^ -9/4 

°3 2 ' 96 ' ' °4 2 ' ' 96 ' ' 

are constants depending on 7. We now fix 7 as 
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4" 4 / 3 -2 + -4- + ,5 + 4 li-T^ >0, (4.5) 



with a < and a E K.. In (|4.5p and the rest of the paper we choose the top sign when we 
are in the Brownian paths setting and the bottom sign when dealing with the two-matrix 
model. It is readily seen that 7 satisfies 

5 - 1673/2 - 1±™ (4-6) 



16 7 3/ 4 2|a|3/2 
and 

7 = 4 ' 4/3 ( lT ^ + 0(<1 " 3> ) (47) 

as a —7- —00 with <r fixed. 

Since 7 depends on the parameters a and <r, the functions and constants Wj, \j,c, Cj 
as well as some functions and constants to be introduced later, depend on a and a 
through 7, although this is not indicated explicitly. Recall that we are interested in the 
behavior of these notions for large negative a. We will add * to the notation to denote 
these functions or constants in the limit a — > —00 with a fixed. 

The properties of the A-functions are listed in the following lemma. 

Lemma 4.2. The functions \j(z), j = 1,2,3,4, defined by (|4.4p have the following 
properties. 

(a) Xj(z) is analytic on IZj, j = 1,2,3,4, and 

Xi,±{x) = A 3 , T (x), x £ (0,oo), (4.8) 

A 2 ,±0c) = A 4 , T (x), x€(-oo,0), (4.9) 

A3,±0*) = a 4,t( z )> z g H c ' °) u (0. ic )- ( 4 - 10 ) 

Hence the function (J J=1 7£j — > C : IZj 3 z 1— >• \j(z) has an analytic continuation 
to a meromorphic function on the Riemann surface 1Z. 
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(b) We have the following symmetry properties 



\j(z) = \j(z) zeKj, j = 1,2,3,4, 

Xi(-z) = X 2 (z), z€C\(-oo,0], 
X 3 (-z) = X 4 (z), zeC\ ((-oo, 0] U [-ic, ic}) . 

(c) As z — > oo within C \ [0, oo), we have 

X 2 (z) = \z^ + (l ± -^j z - *V2 + i _ + 0{z -^ 

A 4 (z) = -§z 3 / 2 + (l ± ^jL^j z + z 1 ' 2 + I + 2DZ" 1 ' 2 + 0{z-\ 

where £ and D are complex constants depending on 7 . 

(d) For z in a neighborhood of the origin we have 

\ 2 ( z ) = G(z)z 2 ^ + H(z)z 4 / 3 + K(z)z 2 , inC\ (-00, 0], 

'L(z), in I U IV, 

< G(z)ujz 2 / 3 + H(z)oj 2 z 4 / 3 + K{z)z 2 , mil, 
G(z)lo 2 z 2 / 3 + H(z)ojz 4 / 3 + K(z)z 2 , in III, 



G(0) = 


32^^ - 16 ^ 2 )' 


G*(0) 


= 0, 


tf(0) = 


3 2 7 > 5 - 16 ^' 


iT*(0) 


= 3 • 4~ 2 / 3 , 


tf(0) = 


1 7 3 / 4 (15 + 16 7 3 / 2 ), 


K*(0) 


1 

" 24' 


L(0) = 


96 7 1 9 / 4 (5 + 48 ^ /2) ' 


L*(0) 


16 

" Y' 



where G(z), H(z), K(z), and L(z) are even analytic functions in a neighborhood 
of the origin with 

(4.11) 

(4.12) 

(4.13) 
(4.14) 

(e) As z — > ic we have that 

X 3 (z) = G{z) + H(z)(z - icf 2 , (4.15) 

where we take the branch cut of (z—ic) 1 ^ 2 along (—ioo,ic) and —n/2 < arg(z-ic) < 
Here G and H are analytic functions in a neighborhood of ic satisfying 

G(ic) = (3 + 16 7 3/2 ), G*(ic*) = 9, (4.16) 

~ 3 1 / 4 ~ 2 

H(ic) = - — (5 + 16 7 3 / 2 )e 3 -/ 4 , H*(ic*) = (4.17) 
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Note that the complete behavior of Xj, j = 1,2,3,4, around the origin, oo, and ±ic 
can be obtained from the symmetry conditions in (b). 

Proof. The proof is straightforward using (|4.4p and Lemma l4.1i □ 

As a — > — oo, the functions Xj converge to A*, i.e., the A-functions associated with 
the parameter 7* = 4~ 4 / 3 . The convergence will be uniform in compact subsets of C. 
For later use, we need the following estimate of the convergence rate. 

Lemma 4.3. There exists a constant g > such that, for every large a, we have 

\Xj(z) - X*(z)\ < £|ar 3/2 max(l,M 3/2 ), z 6 C \ D(0, 5), 

for j = 1,2,3,4, where D(0,6) is an open disk centered at the origin with small radius 
5>0. 

Proof. Due to (|4.7|) . we have 

7 - 7 * = C(|ar 3/2 ), as a -00. 

Together with flJ2]) and this implies 

\Xj(z) - X*(z)\ = 0{\a\~ 3/2 ), as a -00, (4.18) 

where the constant is uniform for z in compact subsets of C. Combining this with 
the behavior of the A-functions at infinity, as given in Lemma 14.21 (c), we obtain the 
lemma. □ 



5 Steepest descent analysis for M 

In this section we perform the steepest descent analysis of the tacnode RH problem ll.il 
As discussed in Sections [23] and E3J we want to study the solution M of this RH problem 
for the following values of parameters 

n = r 2 = 1, 

and 

a 2 . / o~ 



T , HJl ijj-pjj, (5.1) 

as a tends to —00. For the Brownian paths setting it is more convenient to choose the 
plus sign in the formula for t, whereas the minus sign is more appropriate in the two- 
matrix model. The analysis consists of a series of explicit and invertible transformations 

M ^ M (1) ^ M (2) ^ M (3) ^ M (4) ^ M (5) 

of the RH problem. In Section [6] we will use these transformations to prove Theorems 
and 13.31 Since we are interested in the case a — > —00, we may assume that a < 0. 
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5.1 First transformation: M h-> 

This transformation is a rescaling of the RH problem for M. Define 



M^\z-,a)=Ai^(\a\ 1 l 2 M V2 M\~ V2 M\~ V2 w{f^-'^M\( 



Then satisfies the following RH problem. 

Lemma 5.1. The function AT™ defined in (|5.2p /ios £/ie following properties 

(1) M^^z) is analytic for z S C \ Em , where Hm is shown in Figure^ 

(2) AfW /ms the same jump matrices on T,m as M; see Figured 

(3) As z — > oo with z G C \ Sm? aaue 

/\fW(z) = (I + 0{z- 1 )) B{z)A 

x diag ^ e l"l 3 (-^(-2)+^) )e l«l 3 (-^)-^) ;e l«l 3 (^(-^)+^) )e l«l 3 (^)-^)^ ; (5.3) 

where A and B(z) are given in (|1.7p - (|1.8p . and 

(4) AfW is bounded near the origin. 

Proof. This is immediate from RH problem 11.11 and (|5.2p . □ 

5.2 Second and third transformations: i-> M^ 2 ^ i-> 

It is the aim of the second and the third transformation to eliminate the jump matrices 
on Tj, j = 2,3,7,8, see Figure [2j For convenience, we fix the angles ipi = tt/4 and 
<P2 = 7r/3. In a first transformation i— > we erase the lower right block of the 
jump matrices on these rays by moving them to the shifted rays emanating from ic (for 
T2 and T^) and — ic (for and Tg). Thus, we introduce the new rays 

F2,3 = F 2 ,3 + ic, r 7i8 = r 7)8 - ic, 

and define as follows. For k = 2, 3, 7, 8 it is given by 

in the region bounded by T^,, r&, and iM. Here Eij denotes the 4x4 elementary matrix 
with entry 1 at the (i,j)th position and all other entries equal to zero. In the remaining 
regions we simply set = 
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In the next transformation i— )■ we eliminate the remaining upper left block 
of the jumps on Tj, j = 2, 3, 7, 8, by, respectively, moving it to the rays Tj, j = 1, 4, 6, 9. 
This completely eliminates the jump matrices on Tj, j = 2, 3, 7, 8. Thus, for k = 1, 3, 6, 8, 
we define 

for z in the sector bounded by T}. and In the other regions we put 

M (3) _ M (2)_ 

A straightforward check then yields that has jumps on the contour 



E M(3 ) :=1U e™ /4 M U e'^^R U T 2 U T 3 U T 7 U T 8 U [-ic, ic], (5.7) 
as illustrated in Figure [61 and satisfies the following RH problem. 
Lemma 5.2. The function M^ 3 ' as defined in ()5.6p has the properties 

(1) M^\z) is analytic for z G C \ S M ( 3 ) . 

(3) (3) 

(2) M_j_ (z) = M_ (z)J M (3)(z), for z G Tj M ( 3 ), where the jump matrices Jm(3)( z ) are 
shown in Figure [6j Note that we reversed the orientation of some rays. 

(3) As z — > oo, M^ 3 ) has the same asymptotics as ; see (|5.3[) . 

(4) M^ 3 ) is bounded near the origin and near ±ic. 

5.3 Fourth transformation: H> M< 4 ' 

In the fourth transformation we normalize the behavior around infinity using the A- 
functions introduced in Section [4. 2i We define 



M< 4 )(z) = e" |a| ^ (/ - 2iD\a\ 3 E 3>1 + 2iD\afE^ 2 ) M^(z) 

X diag ^ e n 3 ^) )e H 3 ^)^\a\ s X 3 (z)^\a\^\4z)^ ^ ^ 



where i and D are the constants in property (c) of Lemma [ 

Then satisfies the following RH problem with Sj^(4) := £ M (3). 

Lemma 5.3. The function , as defined in (|5.8p . /ias £/ie following properties 

(1) M( 4 )(z) is analytic for z G C \ S M (4) . 

(2) For z G S M (4), we Ziawe M^ 4 (z) = M_ (2) Jjvf( 4 ) where 

J MW (Z) = diag ( e -|a| 3 A 1 ,-(^) je -|a| 3 A 2 ,-W )e -|a| 3 A 3 ,_(,) )e -|a|3A 4 ,_ W ^ 

X J M(3) (z)diag ^l-l 3 Ai, + W )e H 3 A2, + ( 2 ) )e |a|3A3, + (,) )e |a|3A4, + ( 2 )^ _ (59) 
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Figure 6: The jump contour S M ( 3) and jump matrices for M^\z). 



(3) As z — > oo we have 

M^(z) = (I + O (z- 1 )) B{z)A. 

(4) M^ 4 ) is bounded near the origin and near ±ic. 
The jump matrices J M n) (z) are explicitly given by 



(I 0\ 

1 

10 

\0 -1 0/ 



z G (-oo,0), J M (4)(z) 



(5.10) 








1 


o\ 





1 








-1 











\o 








V 



z G (0, oo), 



j mw( z ) = 



/l \ 

10 

1 

-1 el a l 3 ( A4 '+( z )~ A4 '-( z ))/ 



z G (-ic,0) U (0,ic), 
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and 





(*) 


= /-e |a 


3 (Ai(*) 




a| 3 (Ai(2)- 








Jmw 


(*) 


= /-e |a 


3 (A 2 W 


-^Ekfl + e 


a| 3 (A 2 (z)- 


" A4(Z)) ^4,2, 


z e r 4 




Jmw 


(*) 


= / + e |a 


s (A a («) 


-^Eh^ + e 


a| 3 (A 2 (z)- 


" A4W) ^4,2, 






Jmw 


(*) 


= / + e |a 


8 (Ai(*) 


-^E 2A + e 


a|S(Ai(*)- 




zer 9 




Jmw 


(*) 


= / + e |a 


<> 








2 6 f 2 


uf 8 


Jmw 


(*) 


= /-e |a 


3 (A 3 W 


- A4(Z)) S 4 ,3, 






zef 3 


uf 7 



Proof. The explicit formulas of J^f(4) follow from (|5.9p and condition (a) in Lemma 14,21 
To establish the large z behavior of shown in item (3), we first observe from 
condition (c) in Lemma 14.21 and (|5.4p that 

e [a|»(Ai(z)-£(-*)+fc-<) = ! _ 2D| a | 3 (-z)- 1 / 2 + ©(z- 1 ), 
e |o|»(A a (*)-^(*)-fe-<) = ! _ 2j D| a | 3 Z - 1 /2 + 0(z- 1 ), 

e |«|»(As(*)+^(-*)+&-<) = ! + 2D|a| 3 (-zr 1 / 2 + ©(z- 1 ), 
e |o|3(A4(*)+^(*)-fe-<) = ! + 2 £)|a| 3 2- 1 / 2 + ©(z- 1 ), 

as z — > oo. Hence, from the asymptotic behavior of stated in item (3) of Lemma f5.2| 
it is readily seen that 

e -kl^ M (4) (z)diag ^|a|3A l{ ,) )e |a|3A 2 (,) ;e |a|3A3(,) ;e |a|3A 4 ( 2 )^ = ^ + O^ 1 )) B(z)A 

x(l + 2D\a\ 3 diag (-(-z)" 1 / 2 , -z" 1/2 , {-zy 1 ' 2 , z" 1 / 2 ) + dmg{0{z' 1 ))^ , (5.11) 
as z — > oo, where we recall A and -B(z) are defined in (|1.7p - (|1.8p . Note that 

B{z)A (i + 2D\a\ 3 diag (-(-z)~ 1/2 , -z^ 1 / 2 , (-z)" 1 / 2 , z" 1 / 2 ) + diag^Os -1 ))) 

= (J + 2i J D|a| 3 £ 3i i - 2iD|a| 3 £ 4i2 + ©(z^ 1 )) B(a?)A 
This, together with fl52D and (f5TTTD . implies (f5TTUD . □ 

5.4 Estimate of J M (4) on S M (4) 

Let us have a closer look at the jump matrices J M (a) defined in (|5.9|) . The jump matrix 
is constant on (— oo,0) and on (0, oo). On (— ic, 0) U (0,ic) there is a nonconstant (4,4)- 
entry. However, this entry is exponentially small as a — > — oo. Moreover, the decay is 
uniform for z bounded away from the branch points He* , where c* = lim a _ 5 ._ 00 c = 3\/3; 
see (|4.ip and (|4.7p . Also on the other parts of S M (4) the nonzero off-diagonal entries of 
the jump matrices turn out to be exponentially small for large \a\. Again the decay is 
uniform if we exclude small disks around and ±ic* . 
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The preceding statements can be seen as a corollary of the following estimates. We 
denote with D(zq,S) the fixed open disk centered at zq with small radius 5 > and let 
8D(zq,5) stand for its boundary. Recall also that we fixed the angles ip\ = ir/4 and 
(p 2 = vr/3. 

Lemma 5.4 (Estimates for Xj on S M (4)). 

(a) There exist constants cx,C2 > such that 

Re(X 1 (z)-X 2 (z))<-c 1 \z\ 2 / 3 , /or * e(TiUr 9 )\ £>((), <*), (5.12) 

Re(A 1 (z)-A 3 (z))<-c 2 |z| 2 / 3 , for z G (Ti U T 9 ) \ I>(0, 5), (5.13) 

Re(A 2 (z)-Ai(z))<- Cl |z| 2 / 3 , /or z G (r 4 U T 6 ) \ D(0, 5), (5.14) 

Re (A 2 (z) - A 4 (z)) < -c 2 |z| 2/3 , /or z G (r 4 U T 6 ) \ D(0, 5), (5.15) 

/or |a| large enough. 

(b) There exists a constant C3 > suc/i 

Re (X 4 (z) - X 3 (z)) < -c 3 \zf/ 2 , forzef 2 ur 8 \D(ic*,5) 
Re (A 3 (z)-A 4 (z)) < -c 3 N| 3/2 , / rzef 3 uf 7 \D(i C *,^ 

if \a\ is large enough. 

(c) There exists a constant c 4 > such that 

Re (A 4)+ (z) — A 4j _(z)) < — c 4 , 
/or 2 G c* + 5),0) U (0,i(c* — (5)] i/ |a| is sufficiently large. 

Proof. Let us start with claim (a). It is sufficient to prove the first two estimates on Fi, 
since the remaining estimates follow from these ones by the symmetry conditions (b) in 
Lemma 14.21 Moreover, due to Lemma 14.31 and the triangle inequality it suffices to prove 
these estimates for the critical A-functions X*, j = 1,2,3,4. 

Note that T\ = [0, exp i7r/4oo). For large values of z G T\, the estimates follow from 
the asymptotics of A*, which can be obtained by taking a —> —00 in Lemma I4.2f c). 
Using an extra argument one can extend the estimates to Ti \ D(0, 5). This is illustrated 
in Figure El^a). This plots shows the values of ReA*(z), j = 1,2,3,4, for z G Ti. It is 
clearly seen that Re (X*(z) - X*(z)) < 0, j = 2,3, for z G Ti \ D(0,5). This completes 
the proof of (a). 

The proof of (b) is analogous and is based on the plot in Figure E^b) . We recall that 
T 2 = ic* + [0, exp i7r/3oo). Also the proof of (c) is analogous and supported by Figure 
E □ 

Lemma 15.41 has the following immediate corollary. 

Corollary 5.5 (Estimates of J M (4) on S M (4)). 
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Figure 7: Plots of Re\*(z) for j = 1 (solid), j = 2 (dotted), j = 3 (dashed), and j = 4 
(dashdotted). In picture (a), z G Ti, i.e. z = e 7 ™/ 4 ^ with x on the horizontal axis. In 
picture (b), z G I^, i.e. z = ic* + e 7rJ / 3 x with x on the horizontal axis. 

(a) There is a constant c\ > snc/i f/iaf 

J M(4) (*) = / + O (V c iH 3 W 2/3 ) as a -> -oo, 

uniformly for z G S M(4) \(RU [-ic* , ic*] U D(0, 5) U D(±ic* , <5)) . 

(b) There is a constant C2 > suc/i f/iaf 

J M(4 )(z) = I + O (V C2|a|3 ) as a ^ -oo, 
uniformly for z G [i(—c* + 5), i(c* - 6)] \ {0}. 
5.5 Construction of the global parametrix M^> 

If we suppress all entries of the jump matrices for that exponentially decay as 

a — > — oo, we are led to the following RH problem for the global parametrix M^°°\ 

RH problem 5.6 (Global parametrix). We look for a 4 x 4 matrix valued function 
M (°°) that satisfies 

(1) M(°°) is analytic in C \ S M ( m ) , where the contour S M (oo) consists of the real line 
oriented from left to right and the purely imaginary interval [—ic,ic] oriented from 
bottom to top. 
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-2 



Figure 8: Plots of \\ + {z) (solid) and \\_{z) (dashed) for z G (— ic*,ic*). Here z = ix 
with x on the horizontal axis. 

(2) For z G S M (oo), we /taue 



M^Cx) = Mi oo) (x) 



M ( ™\x) = M K f°>{x) 



(oo), 



M^^z) = Mi oo) (z) 
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for x G (0, oo), 



/or x G (— oo, 0), 



(3) As z — > oo, the following asymptotic formula holds 

Af(°°)(z) = (I + 0(z- 1 )) B(z)A. 



(4) We have 



and 



M {co \z) = O (z~ l/Z ) , 
Af(°°)(«) = ((^ic)" 1/4 ) , 



as z — > 0, 



as z — )• ±ic. 



(5.16) 



(5.17) 



/or 2! G (— ic, 0) U (0, ic). (5.18) 



(5.19) 

(5.20) 
(5.21) 



An explicit solution of this RH problem can be built out of the w- functions introduced 
in Lemma l4.1i 
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Figure 9: Image of the map : 7£ i— > C. The solid and dashed lines are the images of 
the cuts in the Riemann surface 1Z under this map. The solid lines represent the branch 
cut of the function k(w) defined in ()5.23p . 



Lemma 5.7. With the functions wj, j = 1,2,3,4, defined in Lemma \4-l\ the global 
parametrix M^°°^ is explicitly given by 



M (oo) ( 



(F^w^z)) F 1 (w 2 (z)) Fx(w 3 (z)) F^w^z)^ 

F 2 (w 1 (z)) F 2 (w 2 (z)) F 2 (w 3 (z)) F 2 {w 4 (z)) 

F 3 ( Wl (z)) F 3 (w 2 (z)) F 3 (w 3 (z)) F 3 (w 4 (z)) 

\F 4 ( Wl (z)) F 4 (w 2 (z)) F±(w 3 (z)) F A {w 4 {z))J 



(5.22) 



where 
Fi{w) 

F 2 (w) 
and 



7 



15/8 



~15/8 



( w + 7 -3/2 )(w _ 7 -3/2 )2) p 3 ( w) 



ij~ 3/8 i -3/2^2 



2y/6K(w) 



w 



-3/8 



/;(/(•) - ( (w 2 - 7 3 )(^ 2 + ^7 3 ) 



1/2 



(^ + 7" 3/2 ) 2 , 



(5.23) 



is defined in the w-plane with branch cut along wi j +(IR + ) U w 3> +((— ic, ic)) U 102,+ (R ) 
(see Figure^) such that 

k(w) = w 2 {l + 0(l/w)), 



as w — ^ +oo. 
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Proof. The jump condition follows from item (a) in Lemma 14.11 and (|5.22p . A straight- 
forward calculation with the aid of item (c) , (d) , and (e) in Lemma 14.11 leads to the 
asymptotic behavior of M^°°\z) in (|5TT9"]) . (15T201) and flOTJ) . □ 

The global parametrix M^°°^ is a good approximation for for values of z bounded 
away from and ±ic*. We will construct local parametrices near these points in the 
next two sections. 

5.6 Construction of the local parametrices around ±ic* 

In the disks D(±ic*,5) around ±ic* we build local parametrices M^ lc ^ respectively. 
The idea is to construct M^ 1 ^ such that it makes exactly the same jumps in the disk 
D(±ic*,5) as M^ 4 ), and satisfies 

M {±ic \z) = (/ + 0(\a\- 3 )) M {oo \z), uniformly for z E dD(^ic*,5) \ S M(4) , (5.24) 

as a —7- — oo. This construction can be done using Airy functions and their derivatives. 
As this construction is standard (cf. |18l I2U] ) and the explicit formulas are irrelevant to 
the proofs of our main theorems, we will omit it here. 

5.7 Construction of the local parametrix around the origin 

In this section we construct a local parametrix in the disk D(0, 5) based on the 

Pearcey parametrix; see RH problem 15.101 below. First, we observe that the (4, 4)- 
entry of the jump for on the imaginary axis is exponentially small, uniformly in a 
neighborhood of the origin; see item (b) of Corollary 15. 5[ Therefore, we may ignore this 
entry in the construction of the local parametrix around the origin which then has to 
satisfy the following conditions. 

RH problem 5.8 (Local parametrix around the origin). We look for a 4 x 4 matrix 
valued function : D(0,5) \ S M (4) — > C satisfying 

(1) M^ ) is analytic on D(0,5) \ S M (4), where S M (4) = S M ( 3 ) is defined in (|5.7p . 

(2) For z £ D(0, 5) n S A /(4), we have 

M i ° ) (z) = M { °\z)J^( Z ), 

with as indicated in Figure [771 

(3) As a — > —oo, the asymptotic formula 

Af(°)(z) = (j + 0(l«r 3/2 )) M(°°\z) (5.25) 

uniformly holds for z 6 dD(0, 5) \ S M (4) . 

(4) M< 0) is bounded near the origin. 
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Figure 10: The jump contour and jump matrices for the local parametrix 
around zero. On the circle the matching condition (|5.25j) is imposed. 



In a first step to solve the RH problem 15.81 we make a transformation that reduces 
the jumps of to constant matrices. Moreover we want the jumps on the imaginary 
axis to disappear. To this end we define 

P( Z ) = M(°)(z)diag ( e -kl 3 AlW 5e -|a| 3 A 2 (,) )e -|a| 3 A 3 (,) ;e -|a|3A 4 (,)^ ^ ^ 

for ±Re z > 0, where C± are two constant matrices given by 

C+ = h, C_ = dia g ri,l,f° 1 J^V (5.27) 

A straightforward check shows that the jumps for M^°\ as indicated in Figure [TTT are 

indeed constant. Furthermore, all jump matrices have the block form ^\ , where 

Q is a non-trivial 3x3 constant matrix. Combined with the asymptotic behavior of 
the A-functions around the origin (see items (b) and (d) in Lemma H~2|) , this leads us to 
propose the following model RH problem. 

RH problem 5.9 (Model RH problem for the local parametrix around the origin). We 

look for a 3 x 3 matrix valued function &(-;p) depending on a parameter p E M that 
satisfies the following conditions. 

(1) is defined and analytic on C\S$, where the contour consists of 6 semi-infinite 
rays as shown in Figure [TM 
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Figure 11: The jump contour and jump matrices for M^- ' around zero. 

(2) For z G £$, we have 

where the jump matrix J$ is constant on each ray and specified in Figure [7 

(3) is z oo and ilmz > 0, we /iaue 



where L± are constant matrices 



-uj 
-1 

-LO 2 





-1 U) 
-1 1 

-1 W 2 



and 0(z; p) is given by 
with 



di&g(-6 2 (z; p), -6 3 (z; p), -6i(z; p)), forlmz > 0, 
diag(-6i(z;p),-d 3 (z;p),-0 2 (z;p)), forlmz < 0, 



%{z- p) = -cj 2k z 4 / 3 + |u, fe z 2 / 3 , fe = 1,2, 3. 



(4) is bounded near the origin. 



(5.28) 



(5.29) 



(5.30) 



At this moment it is not clear where the exact formulation of the asymptotics in item 
(3) comes from. This follows a posteriori from the solution of the RH problem given in 
Lemma 15.111 
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Figure 12: The jump contour and jump matrices for the model RH problem <I>. 



Note that depends on the parameter p through the asymptotics (|5.28p . It turns 
out that we can solve this model RH problem using the Pearcey parametrix that was 
introduced by Bleher and Kuijlaars in the context of random matrices with external 
source [12]. The Pearcey parametrix is the unique solution to the following RH problem. 

RH problem 5.10 (Pearcey parametrix). 

(1) <J> Pc (-;p) is a 3 x 3 matrix valued function depending a parameter p G R, defined 
and analytic on C \ £<j>. 

(2) For z G S^ e , the function <3> Pe has the jump 

where the jump matrix J^ e is constant on each ray and specified in Figure \13l 

(3) As z — > oo and ilmz > ; we have 



d,Pc (z ) = "^V 2 / 8 diag (V 1 / 3 , 1, z 1 / 3 ) L p ± c (/ + 0(z- 2 / 3 )) e 



e(z;p) 



where L^ 6 are constant matrices 





CO 2 














1 








1 














UJ 2 





L Pe = I -1 1 1 | , L^ c = | 1 1 1 | , (5.31) 

and Q(z) given by 

Q( z .p\ = i^s(0i(z;p),e 2 (z;p),e 3 (z;p)), forlmz>0, 
\diag(9 2 (z;p),9 1 (z;p),9 3 (z;p)), forlmz<0, 

and where 9k{z;p) is defined in (|5,30D . 
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Figure 13: The jump contour and jump matrices for the Pearcey parametrix 

(4) <I> Pe (z) is bounded near the origin. 

The above RH problem has a unique solution $ Pe , given in terms of solutions of the 
Pearcey differential equation 

y"'(z) - pp'(z) - zy(z) = 0. 

Since the exact formula of 3> Pe is not relevant for our purposes, we decide not to write 
it down and refer to [121 Section 8.1] for details. 

The next lemma states that RH problem 15.91 can be solved in terms of the Pearcey 
parametrix. 

Lemma 5.11. Let <3? Pe be the Pearcey parametrix as given in RH problem \5. 1U[ the 
solution of RH problem 15.91 is given by 

/0 -1\ 

= <5> Pc (z)- T 1 , (5.33) 
\0 1 / 

where the superscript ~ T stands for inverse transpose. 

Proof. It is straightforward to check that the right-hand side of (|5.33p satisfies RH 
problem 15 .91 Then (|5.33p follows from the uniqueness of the solution to this RH problem. 

□ 

Now, we are ready to construct the local parametrix In view of (|5,26p . we look 

for a parametrix in the form 



X diag ^ e l«l 3 AiW )e |a| 3 A 2 ( Z ) )e |a|3A 3 ( 2 ) 5e |a|3A 4 ( 2 )^ ^ ±Re % > ^ (5^34) 
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where Eq(z) is an analytic prefactor, f(z) is a conformal map near the origin, (fi(z) is an 
analytic function and C± are the constant matrices shown in ()5.27p . The parameter p 
will also be chosen to be dependent on z. 

Recall the functions G(z), H(z), and K(z) introduced in item (d) of Lemma [4.21 We 
first define 

3/4 



f(z) = f(z;a) = \a\^ (^H(z)j 
By (|4.12p . it is readily seen that 

'4 



z, 



z G D(Q,6). 



1 9/4 



3/4 



H(0) z + 0(z i 



(5.35) 



as z — > 0. Hence, / is a conformal map in a neighborhood of the origin. We then modify 
the contours (S M ( 4 ) \ iU.) n D(0, 5) if necessary in such a way that / maps them into 
£$. Next we set 



p{z) = p(z; a) = 2\a\ 



>G(z)z 



2/3 



and 



(z: a 



/(z) 2/3 

-|a| 3 A 4 ( Z ) ; 
-\a\ 3 M(z) 



z g D(o,<y), 



Rez > 0, 
Rez < 0. 



(5.36) 



Also the function 4>(z) is analytic on D(0,5), which follows from (|4.10p . With these 
definitions of f(z) and p(z), it is easily seen from items (b) and (d) in Lemma 14.21 that 



h(f(z); P (z)) 



>| 3 (A 3 




- K{z)z 2 ) , 


for 


z e D(p, 5) n I, 


< \a 3 (A 4 




- K{z)z 2 ) , 


for 


z g D(p,8)nii, 


Ja| 3 (Ax 




- i^Wz 2 ) , 


for 


z g D(Q,S) n (in u IV) 


>| 3 ( A i 




- K{z)z 2 ) , 


for 


z g D(o,5) n (lull), 


< a 3 (A 4 


» 


- K{z)z 2 ) , 


for 


z g D(o,5)r\in, 


>| 3 (A 3 




- i^Wz 2 ) , 


for 


z G D(o,5)niv, 


a| 3 (A a (* 


)- 


K(z)z 2 ) , 


for 


z G D(0,S) \ (-5,0}. 



(5.37) 



Finally, the matching condition (|5.25p leads us to the definition of the prefactor Eq(z) 
as 



Eo(z) = M^(z)(C ± )- 1 

-iV67rePW a /8-|ol 8 ^W* 2 (L±)-i diag (/(z)- 1 ^ i, f(z) 1 / 3 ) 

1 



±Rez > 0, 
(5.38) 



where L± is defined in ()5.29p . We then have the following lemma. 
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Lemma 5.12. Eq(z) has an analytic continuation to D(0,5). 



Proof. Eq(z) is clearly analytic in D(0, 6) \(1U iR). For z G iR \ {0} we have that 
Eq^{z) = Eq^(z) as a consequence of (|5. 18[) and ([5.270 . There is no jump on R + 
either, because of (|5,16p and the observation 




With a bit more effort one can also check that the jump on R~ is trivial. Hence Eq(z) 
is analytic in the punctured disk D(0, 5) \ {0}. It follows from the behavior of the global 
parametrix M(°°) around zero that the singularity of Eq(z) at the origin is removable. □ 

It is now straightforward to check that the jumps of on D(0,5) n S M (4) and 

its large \a\ behavior on dD(0,5) (with the aid of f)5.28|) and ([5.370 ) are indeed of our 
required form. In summary, we have the following lemma. 

Lemma 5.13. The matrix valued function M^ ' defined in ([5,340 satisfies conditions 
(l)-(4) of RH problem \5M 



5.8 Final transformation: M< 4 ' ^ M^> 

Using the global parametrix Af' 00 ' and the local parametrices M( ±ic ) and Af(°) , we define 
the fifth transformation \— > as follows 



' ' mW(z){M(°°\z)) 1 for z G C\ (S M( 4) UD(0,8) UD{±ic*,5)), 



M&{z) {M^(z))- 1 for z G D(0,S) \ S M(4) , 

k M( 4 )(z) {M^ ±lc \z)y l for z G D(±ic*,5) \ S A/(4) . 

(5.39) 

Then is defined and analytic outside of S A/ (4) and the three disks around and 
±ic* , with an analytic continuation across those parts of E M (4) where the jumps of the 
parametrices coincide with those of M^. What remains are the jumps on a contour 
S M ( 5 ) that consists of the three circles around 0, ±ic* , the parts of Ti, I^, T3, T^, Tg, 
Ft, Tg and Tg outside of the disks, and (i(—c* + 5),i(c* — 5)). The circles are oriented 
clockwise. Then satisfies the following RH problem. 

RH problem 5.14. 

(1) M^ 5 ) is defined and analytic in C \ £ M (s). 



(5) 

(2) For z G Ew ; we have 



M®(z)=M ( ?\z)J M(B) (z), 
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where 




= < 



i 



for z£ dD(0,6), 
forze dD(±ic*,5), 
for z G (— i<5, iS), 
elsewhere on E M (5). 



(3) As z — > oo, we have 



M®(z) =I + 0{l/z) 



Furthermore, the jump matrix J M (S) tends to the identity matrix on S M ( 5 ) as a — > 
— oo, both uniformly and in L 2 -sense. Indeed, from the matching conditions (|5.24p and 
(|5.25p . we see that 



as o 4 — oo, uniformly on the circles around and He*. On the remaining parts of 
^m( 5 )> ^m( 5 ) ^ s uniformly exponentially small (see Corollary I5.5|) . Then, as in [332 [20], 
we conclude that 



as a — > — oo, uniformly for z in the complex plane outside of S M ( 5 ). 

The estimate (|5.4ip is the main outcome of our Deift-Zhou steepest descent analysis 
for RH problem 1 1.11 We will use it to prove our main results in the next section. 

6 Proofs of the main theorems 
6.1 RH formula for the Pearcey kernel 

Before coming to the proofs of the main theorems, we state the following lemma giving 
an expression of the Pearcey kernel in terms of the Pearcey parametrix. 



J M(5) (z) = I + 0(\a\-^ 2 ) 



(5.40) 




(5.41) 
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Lemma 6.1. Let K Pe (x,y; p) be the Pearcey kernel defined in 

(I 



we have 



K Pe (x,y;p)= { 



27T2 (^X ~ 


HI 


l 




2ni(x- 


-y) 


1 




2ni(x- 


-v) 


1 




2ni(x- 


-y) 



-1 1 o)<S>l c (y;p)- l <Z>l c (x;p) 



1 1 )^(y;p)- 1 ^l e (x;p) 



1 1 0)<!> p + c (y;p)- 1 ^ + e (x;p) 



-1 l) &?(y; py^^ix; p) , for x, y < 



1 | , for x, y > 0, 
1 | , /or x > 0, y < 0, 

u 

I , /or x < 0, y > 0, 

u 

1 N 



where <^ Pe (-;p) is t/ie i/je unique solution of RH problem [h.l(A 
Proof. See [121 Section 10.2]. 

6.2 Proof of Theorem I2HEI 



□ 



First we prove ()2.10p . We will focus on the case u, v > 0, or equivalently, x, y > 0. Other 
cases can be proved similarly. We start from (|2.6p . and the strategy is to express this 
kernel in terms of instead of M by unfolding all transformations M i— > i— > 

( 2 ) i-> M^ 3 ) h> M^ 4 ) i-> M( 5 ) of the steepest descent analysis. 
From the first transformation M \— > in (15.21). it follows 



o?K t&c [a 2 u,a 2 v;-±a 2 ,\a\ ( 1 + 



a 



2|a| 3 / 2 



2-ni{u — v) 



1 1 0)M^(v;a) ^Mj lj (it; a) 



(i) 




1 

w 



where, since we work in Brownian paths model, we chose the top sign in ±o~. The 
transformations M^' \— > M^ 2 ) i— >■ M^ 3 ^ in (|5.5p and ()5.6p leave this formula essentially 
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unaffected. Applying the transformation i— > in (|5,8|) . however, yields 



a 2 K tac a 2 u,a 2 v;a,-±a z ,\a\ I 1 + 



i „2 



2|o| 3 / 2 



2m(u — v) 



/ e -|a| 3 Ai, + («)\ 



-| a l 3 <^3,+0) 

V o / 



In the next step we unfold the transformation i-> M^ 5 ) in (|5,39p . Assuming that 
< tt, v < 5, it follows from ()5.34p that 



a 2^tac ( a 2 n;a 2 u; _l a 2 ) | a | ( 1 + 



(7 



2|a| 3 / 2 // 27ri(u-u) 



x(-l 1 o)^^ 1 J )e ( v )-^(v^ 



V 



/1\ 



1 

w 



Now we fix x, y > and take 

u = 2- 1 / 2 \a\- g / A x, 



2- 1 /2| a |-9/4 2 



so that < u, v < 5 for \a\ sufficiently large. Under this change of variables, it follows 
from QOSD and (fl~T2l that 

f(u) -4 x, /(«) -> y, 
as a -4 -oo. Also (|536l) . (IQ5|) . (I4TT21) . (jillT) . and imply 

p(lt) -4 cr, p(v) -)■ (7, 
as a —7- — oo. Furthermore, by standard considerations it follows that 



Af ® (v; ay 1 M^ {u;a)=I + 



v — u 



l| 3 / 2 



I + O 



l| 15 /4 J 



as a — > — oo, uniformly for x and y in a compact subset of M. Observe also that 
Eq(u) = 0(|a| 3 / 4 ) as a — > — oo; see (|5.38|) and (|5.2U|) . The same bound holds for 
Eq(u) , Eq(v) and Eq(v) , so that we find 



E (v)' l E (u) = I + O (\u - v\\a 



I + O 



\x - y\ 
|a| 3 / 4 
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Combining all these results, gives 

l im ——}—-^ K tac ( —pJ——, —=^——\ —ha 2 , \a\ ( 1 + G 



^-oo y/2\a\V 4 Vv / 2|a| 1 /4' V2|a|V4' 2 >ii^ ' 2 |a| 3 / 2 

1 



2iri(x — y) 

An appeal to (|5.33|) yields 



1 1) ^(y^y^ix-^) [ 



lim ( ^ 77 , ; -\a\ \a\ ( 1 + ° 



^-00^10)1/4 VV2H 1 / 4 ' v^lal 1 / 4 ' 2 V 2|a|3/2 

' T 1 O^foaf^fca)- 2, I -1 



2iri(x — y) 







Finally, by taking the transpose on both sides of the above formula, (|2.10p follows from 
for the case x, y > 0. 

Next we prove (|2.1ip from symmetry considerations. We start with the following 
symmetry relation 

M- T «;s,*) = (_° 2 

for s,t£R, which can be checked from RH problem II .li Here, I2 denotes 2x2 identity 
matrix. Using (|2.5p this leads to the following symmetry property of the tacnode kernel 

K^ c (u,v;s,-t) = K tac (v,u;s,t). 

Given this, (|2.1ip is immediate from ()2.10l) . 
This completes the proof of Theorem 12. 31 

6.3 Proof of Theorem [3731 

We start with a lemma that establishes a symmetry property of the Pearcey parametrix. 
Lemma 6.2. Let $ Pe (-; p) be the unique solution of RH problem \5.1(A Then 

diag (— i, 1, i) & Pe (iz; —p)Bj = <& Pe (z; p), for z in the j-th quadrant, (6-2) 

where 

-l\ / 1 0\ /0 1 0\ /0 -1 N 

Bi = I -1 , B 2 = -1 , £3=|00l], £4=[01 

010/ V 01/ \1 0/ \1 
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Proof. We have to check that the left-hand side of (|6.2|) satisfies RH problem 15. 1UI It is 
straightforward to check the jump conditions in item (2). Checking the asymptotics is 
a bit more cumbersome. To that end the following observations are useful 



h(iz; -p) 
h(iz; -p) 
%{iz; -p) 





P), 


for z G / U /// U IV, 


i(*; 


p), 


for z G II, 


i(z; 


P), 


for z G / U III U IV, 




P), 


for z G //, 


2(z; 


P), 


for z G / U III U IV, 


s(z; 


b), 


for z G II, 



and also 

e ±«/6(^)Tl/3 = z Tl/3 ; 
Ti(^) Tl/3 = ^ TV3 , 

We omit the details here. 



for zelUHIU IV, 
for z G //. 



□ 



The proof is similar to that of Theorem 12.31 and we will focus on the points where 
both proofs differ and leave some details to the reader. Note that we are now working 
in the context of the two-matrix model, so we choose the bottom sign in ±cr; cf. (|4.5p 
and ([5ZED . 

Following the same ideas as in the proof of Theorem 12.3} we obtain 



a 2 K" 



2 2 

a u,a v; 



i „2 



a , a II 1 



1 



(-1 1 0) 



2|a| 3 / 2 y ) 2-Ki(u - v) 



x Afi 5) (w)£b(*«) 



o 



o 

4>(iv) 



1 




for < \u\, \v\ < 5. 

Now we fix x, y ^ and take 

u = 2~ 1 / 2 \a\- 9 / 4 x, v 

so that \u\, \v\ < 5 for \a\ sufficiently large. Hence 

f(iu) -)• zx, /(it; 

as a — > — oo and 

p(iu) -> -cr, 



2- 1 /2| a |-9/4 2/; 



-cr, 
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as a —7- — oo. Furthermore, 



M| 5) (iu; a)- l Mf (iv; a) = I + O 1^—^ } . 
as a — > — oo, uniformly for x and y in a compact subset of R. Also 



^o(^)- 1 J B (w) =/ + of|u-u||a| 3/2N ) =/ + C 



|g - y\ 

|a| 3 / 4 

These results, together with (|5,33p and taking the transpose, imply 
lim -Jl^jcxI J.., --U\\a\(l 



«->-oo y/2\a\V 4 V^la) 1 ^' ^2| a |i/4' 2 2 | a | 3 / 2 
Finally, we obtain Theorem 13.31 by applying Lemma 16.21 and (|6. lj) to the above formula. 
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